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This paper presents a critical review of particle production in an uniform electric field and 
Schwarzchild-like spacetimes. Both problems can be reduced to solving an effective one-dimensional 
Schrodinger equation with a potential barrier. In the electric field case, the potential is that of an 
inverted oscillator —x 2 while in the case of Schwarchild-like spacetimes, the potential is of the form 
— 1/x 2 near the horizon. The transmission and reflection coefficients can easily be obtained for both 
potentials. To describe particle production, these coefficients have to be suitably interpreted. In 
the case of the electric field, the standard Bogoliubov coefficients can be identified and the standard 
gauge invariant result is recovered. However, for Schwarzchild-like spacetimes, such a tunnelling 
interpretation appears to be invalid. The Bogoliubov coefficients cannot be determined by using an 
identification process similar to that invoked in the case of the electric field. The reason for such a 
qq , discrepancy appears to be that, in the tunnelling method, the effective potential near the horizon 

■ is singular and symmetric. We also provide a new and simple semi-classical method of obtaining 

Q\ ' Hawking's result in the (t, r) co-ordinate system of the usual standard Schwarzchild metric. We give 

a prescription whereby the singularity at the horizon can be regularised with Hawking's result being 
recovered. This regularisation prescription contains a fundamental asymmetry that renders both 
sides of the horizon dissimilar. Finally, we attempt to interpret particle production by the electric 
field as a tunnelling process between the two sectors of the Rindler metric. 
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I. INTRODUCTION AND SUMMARY 

o 

1 In this paper, we present a critical review of particle production in Schwarzchild-like spacetimes and in an uniform 
electric field in Minkowski spacetime. Our emphasis here is on the tunnelling picture used to interpret results obtained 
OO , by other (presumably more reliable!) methods. 

Consider the problem of a scalar field propagating in flat spacetime in an uniform electric field background. The total 
particle production due to the presence of the electric field upto the one-loop approximation is correctly calculated 
by the gauge invariant method proposed by Schwinger |Q. The same problem can be reduced, in a time dependent 
1 I gauge, to an equivalent Schrodinger equation with an inverted harmonic oscillator potential. Such an equation can 
be solved exactly using the standard flat spacetime quantum field theoretic methods. Since the problem is explicitly 
time dependent, the vacuum state at t — ■> — 00 and at t — > 00 are not the same. The Bogoliubov coefficients between 
the "in" and "out" vacua are easily calculated and the total particle production turns out to be the same as calulated 
by the Schwinger method. However, if a space dependent gauge is used to describe the same field, the vacuum state 
of the field remains the same for all time and hence no particle production can take place. To recover the standard 
result, the tunnelling interpretation is introduced. This interpretation is useful since it provides a dynamical picture 
of particle production (see, for example 0) and is the only way in which the standard gauge invariant result can be 
recovered in a time independent gauge. In this paper, we attempt a tunnelling description for both time dependent 
and time independent gauges. The method of complex paths, enunciated by Landau in |l0| ], is used to calculate the 
transmission and reflection coefficients for the equivalent quantum mechanical problem. Then, an interpretation of 
these coefficients, in order to obtain the standard gauge invariant result is provided. 

Let us next consider Schwarzchild-like spacetimes like the usual black hole spacetime, the Rindler spactime and 
the de-Sitter spacetime. In the standard black-hole spacetime, particle production was obtained by Hartle and 
Hawking Q using semi-classical analysis. In this method, the semi-classical propagator for a scalar field propagating 
in the Schwarzchild spacetime is analytically continued in the time variable t to complex values. This continuation 
gives the result that the probability of emission of particles from the past horizon is not the same as the probability 
of absorption into the future horizon. The ratio between these probabilities is of the form 
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P [emission] = P [absorption] e ^ E (1) 

where E is the energy of the particles and (3 — 1/8%M is the standard Hawking temperature. The above relation is 
interpreted to be equivalent to a thermal distribution of particles in analogy with that observed in an atom interacting 
with black body radiation. In the latter case, the probability of emission of radiation by the atom is related to the 
probability of absorption by the atom by a similar relation as given above. Now, in Hawking's derivation, the 
Kruskal extension is of vital importance in obtaining the thermal spectrum. In this paper, we propose an alternate 
derivation of Hawking radiation without using the Kruskal extension. We consider the semi-classical propagator 
in the standard Schwarzchild co-ordinates near the horizon and then use a prescription to bypass the singularity 
there. This prescription gives the same result as that obtained by Hawking and can be used in all spacetimes with 
a Schwarzchild-like metric. We then attempt to create a tunnelling picture on the lines of that done for the electric 
field. The equivalent Schrodinger potential in this case is of the form — 1/x 2 near the horizon. The transmission and 
reflection coefficients can easily be obtained using Landau's method of complex paths. To obtain particle production, 
these coefficients have to be suitably interpreted. In the case of the electric field, the standard Bogoliubov coefficients 
can be identified by recasting the normalisation rule between the transmission and reflection coefficients, namely 
T + R = 1, in such a way that the standard gauge invariant effective lagrangian result is recovered. However, in 
the case of the Schwarzchild-like spacetimes, such a tunnelling interpretation appears to be invalid. The Bogoliubov 
coefficients cannot be identified using the transmission and reflection coefficients in the same way as was done in 
the case of the electric field since the spectrum of created particles is then not thermal as Hawking's semiclassical 
method suggests. But reinterpreting these coefficients also does not give the standard result. The reason for such a 
discrepancy appears to be that in the tunnelling method, the effective potential near the horizon is symmetric and 
singular and hence both sides of the horizon appear to be the same. 

Finally, we attempt to link particle production by an uniform electric field with processes occuring in the Rindlcr 
frame. We propose that tunnelling between the two Rindler sectors gives rise to the production of real particles. We 
do this entirely in a heuristic manner and show that this tunnelling process between th e tw o Rindler sectors gives 
rise to the exponential factors in the expression for the effective lagrangian (see Eqn. ( |1 04| ) . But the coefficients 
multiplying the exponential factors, however, must be determined by solving the problem exactly. 

The layout of the paper is as follows. In section (|n|), we will discuss first the calculation of the reflection and 
transmission coefficients using Landau's powerful method of complex paths for the potentials — x 2 and —1/x 2 . Then, 



in section (III) we apply the results of the previous section to particle production in an electric field. In section (IV), 
the semi-classical derivation of Hawking radiation without taking recourse to the Kruskal extension will be extensively 
discussed. Next, the reduction of the problem to an effective Schrodinger equation and tunnelling at the horizon will 
be considered. Finally, in section (^), we attempt to link particle production in the electric field to processes occuring 
in the Rindler frame. 



II. FACETS OF TUNNELING 



In this section wc briefly review the basic concepts of semi-classical quantum mechanics in one dimension and 
formally describe the tunnelling process. We then apply the formalism to two potentials, namely, Vi(x) = —x 2 and 
Vi{x) = —1/x 2 , and calculate the transmission and reflection coefficients for both. 



A. Semi-classical limit of Quantum Mechanics 

Consider a simple one dimensional quantum mechanical system with an arbitrary potential V{x) where x denotes 
the space variable (see Ref. [[l0| for details). To describe the transition of the system from one state to another, 
we first solve the corresponding classical equations of motion and determine the path of transition. This path is, in 
general, complex since many processes like tunneling through a potential barrier cannot occur classically. Therefore, 
the transition point qo where the system formally makes the transition is a complex number determined by the classical 
conservation laws. Then, the action S for the transition from some initial state Xi to a final state Xf given by 

S(x f ,Xi) = S(x f ,q ) + S(q ,Xi) (2) 

is calculated. Here, S(qo,Xi) is the action for the system to move from the initial state Xi to the transition point qo 
while S(xf,qo) is that to move from go to Xf. The probability P for the transition to occur is given by the formula 

P~ exp (~ha[S(xf,q a ) +S(q ,Xi)] \ . (3) 
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The above formula is valid only when the exponent of the exponential is large. Further, if the potential energy has 
singular points, these must also be considered as possible values for qo- If the position of the transition point is not 
unique, then it must be chosen so that the exponent in Eqn. (^) has the smallest absolute value but still must be 
large enough so that the above formula be valid. 

If the transition point qq is real, but lies in the classically inaccessible region, then the above formula gives the 
transmission coefficient for penetration through a potential barrier, while if the transition point is complex, it solves 
the problem finding the over the barrier reflection coefficient. The ~ sign in the above formula is used since the 
coefficient in front of the exponential is not determined. This can be determined by calculating the exact semi- 
classical wave functions. Generally, it is desirable to find the ratios of two different transitions so that this coefficient 
does not matter. 

The physics of the tunnelling and the "over the barrier" reflection processes are very different. In the tunnelling 
process, the semi-classical analysis gives a transmission coefficient that is an exponentially small quantity with the 
corresponding reflection coefficient being unity. In contrast, in the "over the barrier" reflection process, just the reverse 
is obtained. The transmission coefficient is unity while the reflection coefficient is an exponentially small quantity. 
Both these processes will be encountered when the electric field is studied in different gauges. 

We will now review the method of calculating the transmission and reflection coefficients for a typical quantum 
mechanical problem using the method of complex paths for a general potential V(x). 



1. Description of the method of complex paths 

Consider the motion of a particle of mass m in a region characterised by the presence of a potential V(x) in one 
space dimension. The problem is to calculate the transmission and reflection coefficients between two asymptotic 
regions labelled L and R where the semi-classical approximation to the exact wave function is valid. After identifying 
these regions and writing down the semi-classical wave functions, definite boundary conditions are imposed. The usual 
boundary conditions considered are that in one region, say L, the wavefunction is a superposition of an incident wave 
and a reflected wave while in the second region R, the wave function is just a transmitted wave. Then, a complex path 
(in the plane of the now complex variable x) is identified from R to L such that (a) all along the path the semiclassical 
ansatz is valid and (b) the reflected wave is exponentially greater than the incident wave at least in the latter part of 
the path near the region L. The transmitted wave is then moved along the path to obtain the reflected wave and thus 
the amplitude of reflection is identified in terms of the transmission amplitude. Having done this, the normalisation 
condition is used i. e. the sum of the modulus square of the transmission and reflection amplitudes should equal unity, 
to determine the exact values of the transmission and reflection coefficients. 

For a given potential, the turning points qo (or transition points) are given by solving the equation 



p(q ) = y/2m(E-V(q )) = (4) 

where p{x) is the classical momentum of the particle and E is the energy of the particle. In general, go is complex. 
At these points, the semi-classical ansatz is not valid since the momentum is zero. Further, the potential can possess 
singularities. At these points too, the semi-classical approximation is invalid. Therefore the contour between the two 
regions should be chosen to be far away from such points. In general the contour will enclose them. Therefore, the 
relation between the transmission and reflection amplitudes is determined by taking into account the turning points 
and the singularities of the potential. 

The Schrodinger equation to determine the wave function tp of the particle is 

_tfd^ = (E _ v{x)) ^ (5) 
2m ax 

Referring to flo|| , the semi-classical wave function, in the classically allowed region where E > V(x), is given by the 
formula 

yj = dp- 1 ' 2 exp ^ p{x)dx^j + dp' 1 ' 2 exp (- l - J p{x)dx^j (6) 

while in the classically inaccessible regions of space where E < V(x), the function p{x) is purely imaginary and the 
wave function is now given by the relation 



ij = dM-^cxp (~ J \p{x)\dx\ + CM- 1 ' 2 exp Q J |. 



\p(x)\dx). (7) 
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The condition on the potential for semi-classicality of the wave function to be valid is 



d 
dx 



p(x) 



or, in another form, 



nk\F\ 

1pY~ 



«1, 



« 1, 



dx 



(8) 



(9) 



It ought to be noted that the accuracy of the semi-classical approximation is not such as to allow the superposition of 
exponentially small terms over exponentially large ones. Therefore, it is inapplicable in general to retain both terms in 
Eqn. (0). We will consider a few cases of interest in this paper and refer the reader to for a exhaustive discussion 
along with suitable illustrative examples. 

Consider the case in which the semi-classical condition (^) holds in the regions x — » ±oo. As x — > — oo, the wave 
function is assumed to be a superposition of incident and reflected waves and is written in the form, 



1p =p 



exp (~ j p(x)dx^j +C 2 p 1/2 exp ^--^ Jp 



(x)dx 



where the incident wave has unit amplitude while the reflected wave has amplitude given by Ci- As i 
wave function is assumed to be a right moving travelling wave, 



(10) 



-oo, the 



c 3 I i 

— exp 

VP 



pdx 



(11) 



The method of complex paths can now be applied on the function (pd|). The contour is chosen either in the upper or 
lower complex plane such that the reflected wave is always exponentially greater than the incident wave along that 
part of the path near the region x — > — oo. If this is satisfied along one of the contours then C2 is determined in terms 
of C3. To carry out the above procedure however, the exact semi-classical wave functions as x — > ±00 have to be 
determined. This will be done explicitly later for the relevant cases. 

A different case arises when the semi-classical ansatz holds in the vicinity of the origin x — rather than at x = ±00. 
The boundary conditions assumed in this case are the same as above with the condition x — > 00 replaced by x — > + 
and x — *■ —00 by x —> 0~ . Here, the required contour is about the origin and is chosen to be small. But it must still 
be large enough so that the reflected wave is much larger than the incident wave along the latter part of the contour 
near the region x < 0. 

In the above cases the method of complex paths gives the exact transmission and reflection amplitudes. But, 
in certain cases it is enough to assume that the transmission and incident amplitudes are equal to unity while the 
reflection amplitude is exponentialy damped and consequently very small. Here, the "over the barrier" reflection 
coefficient for energies large enough so that the reflection coefficient is exponentially small, has to be determined. 
In this case, the condition E > V{x) is always satisfied. Therefore, the transition point go at which the particle 
reverses its direction is the complex root of the equation V(qo) — E. Let go lie in the upper complex plane for 
definiteness. Now, the amplitudes of the incident and transmitted waves are equal (both are set to unity within 
exponential accuracy). To calculate the reflection coefficient, the relation between the wave functions far to the right 
of the barrier and far to the left of the barrier must be determined. The transmitted wave can be written in the form 



1 / 1 
exp 



VP 



pdx 



(12) 



where X\ is any point on the real axis. We follow the variation of tpT along a path C in the upper complex plane which 
encloses the turning point qo. The latter part of this path must lie far enough to the left of go so that the error in the 
semi-classical incident wave is less than the required small reflected wave. Passage around qo only causes a change in 
the sign of the root \J E — V(x) and after returning to the real axis, the function ipT becomes the reflected wave ipR- 
Going around a complex path in the lower complex plane converts ipT into the incident wave. Since the amplitudes 
of the incident and transmitted waves may be regarded as equal, the required reflection coefficient is given by 



R = 



exp 



Im / pdx 



C 



(13) 



Now we can deform the contour in any way provided it still encloses the point go- I n particular, the contour can be 
deformed to go from x\ to go and back. This gives 
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R = exp 



-Im 



pdx 



(14) 



Since p(x) is real everywhere, the choice of x\ on the real axis is immaterial. The above formula determines the above 
the barrier reflection coefficient. It must be emphasised that to apply the above formula the exponent must be large 
so that 1 — R is very nearly equal to unity. 

Finally consider a situation where the amplitudes of the reflection and incident wave are equal. The transmission 
coefficient is now an exponentially small quantity. This case corresponds to the standard tunnelling process. The 
problem is characterised by the presence of real turning points between which lies the classically forbidden region 
where the energy E < V(x). For definiteness, let there be two real turning points labelled g_ and q + . The potential 
V(x), in the immediate vicinity of the turning points of the barrier, is assumed to be of the form 



E-V(x)^F (x-q ± ), 



F = - 



dV 
dx 



(15) 



This assumption is equivalent to saying that the particle, near the turning points, moves in a homogeneous field. 
With this assumption, the amplitude of transmission C3 is given by (refer to Ref. |10| page 181) 



C 3 



exp 



\p(x)\dx 



(16) 



The transmission coefficient is then given by | C3 1 2 . The above formula holds only when the exponent is large. In 
the derivation above, we have assumed that the semi-classical condition holds across the entire barrier except in the 
immediate vicinity of the turning points. In general, however, the semi-classical condition need not hold over the 
entire extent of the barrier. The potential, for example, could drop steeply enough so that Eqn. ( |T5| ) is not valid. In 
these cases, the exact semi-classical equations have to be determined before applying the method of complex paths. 
The cases encountered in this paper all satisfy Eqn. (p^|). 

We now apply the above results to two potentials. The first is the well known inverted harmonic oscillator potential 



Vl(x) 



-g\x with <7i > while the other is \>2(x) = —gifx with 52 > 0. The first potential arises when the 



propagation of a scalar field in a constant electric field background is studied. The second potential arises when the 
propagation of a scalar field in Schwarzchild-like spacetimes is considered in the vicinity of the horizon. 



2. Application to the potential Vi(x) 



Consider the potential given by 



Vi(x) 



-gix 



(17) 



where g\ > is a constant. This potential is the inverted harmonic oscillator potential and is discussed extensively in 
many places (see for example || , |[o| ) . We will follow the semi-classical treatment given in Ref. jl0| and briefly review 
the calculation of the reflection and transmission coefficients for both the tunneling and over the barrier reflection 
cases. 

The semi-classicality condition (||) for the above potential is, 



hgix 



gix 2 ] 3 / 2 



< 1 



(18) 



where m is the mass of the particle and E\ is its energy. The above condition definitely holds for large enough \x\ 
and for any value of E%, either positive or negative. Therefore the motion of a particle moving under such a potential 
is semi-classical for large enough \x\ and hence holds as x — > ±00. 

Since the motion is semi-classical for large \x\, we can expand the momentum p(x) as, 



p(x) = \Jlm (Ei + gix 2 ) w \f 2mg 1 ^x + 



2gxx 



(19) 



Using Eqn. (|l9|), the semi-classical wave functions can be written as follows with the following boundary conditions. 
As x — > 00, we assume that the wave function is a right moving travelling wave ipR while as x — > —00, it is a 
superposition of an incident wave of unit amplitude and a reflected wave given by tpL- Therefore, we have, 
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4>R = C 3 ? 



where we have made the definitions 



-00) 
-oo) 



(^) 



1/4 



1 



£1 




Si 



(20) 
(21) 

(22) 



Following the variation of Eqn. (|20|) around a semi-circle of large radius p in the upper half plane of the now complex 
variable £, we obtain 



C 2 = -iC 3 exp (— 7T£i) 
The conservation of particles is expressed by the condition that, 

\C 3 \ 2 + \C 2 \ 2 = 1 
the transmission coefficient is, 
1 



From Eqn. (E3h and Eqn. (E 



T=\C 3 \ 2 = 



1 



3 — 27T£i 



1 + e -k*V 2m /9iEi 



while the reflection coefficient is 



R 



\C 2 \ 2 



1 + e ~^ 7r y /2m/9lEl 



(23) 
(24) 

(25) 
(26) 



Note that the passage through the lower half complex plane to determine C 2 is unsuitable since on the part of the 
path — 7r < 4> < — 7r/2 , where (f> is the phase of the complex variable £, the incident wave (first term in Eqn. (pl|)) is 
exponentially large compared with the reflected wave. The above formula holds for all energies E\. This is because, 
even for negative energies, the semi-classical wave functions given in Eqns. (|2(], ^l|) are exactly the same with the 
boundary conditions being fully satisfied. 

If E\ is large and negative, Eqn. (|25|) gives T w e - 7r \ /2,n /9 1 \ El \/ h and thus R ~ 1. This is in accordance with the 
formula in Eqn. (|ll]). To apply Eqn. (p_6|) it is necessary to calculate the turning points first. The real turning points 
are qo = —^J\Ei\j g\ and qi — y/\Ei\/g±. Therefore, 



exp 



\p(x)\dx 



exp ( ~—TT V / 2m/g 1 \E 1 \ 



\Ei\ 
9i 



dx 



(27) 



This gives the same answer. 

We can calculate the over the barrier reflection coefficient using Eqn. ( |l4| ) for E\ large and positive. The turning 
points now are given by solving the equation p(qo) — with the condition that E\ > V\{x) always. Since E\ > 0, the 
turning points are go = ±*\/-Ei/<7i- Choosing the positive sign for go and setting x\ = 0, the integral in Eqn. ( |l~4] ) is 
evaluated as follows: 



p(x)dx = y/ 2mg[ / \J E\jg\ + x 2 
Jo 

= iy/2m/g\Ei I \J\-y 1 = -in v / 2m/g 1 E 1 
Jo 4 



Therefore, 



R = exp ( — — -K\j2ml g\E\ 
\ n 



(28) 



(29) 



The above formula can also be obtained directly from Eqn. ( p6| ) by neglecting the exponential term compared to unity 
which means that the energy has to be large. 
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3. Application to the potential Vz(x) = —gijx 1 



Consider the potential given by 



Vx(x) = -f 2 (30) 



where g 2 is a positive constant. The potential has a singularity at the origin. This potential arises when the effective 
Schrodinger equation is calculated for Schwarzchild-like spacetimes. This aspect will be dealt with in later sections. 
The semi-classical condition (Q) for this potential takes the form 



hg 2 1 



< 1 (31) 



where E 2 is the energy. It is clear that the above relation holds for large \x\. It also holds for small |a;| if \J2mg 2 S> h. 
Notice that the quasi-classicality condition for small \x\ is independent of the sign and magnitude of the energy E 2 . 
For this potential, we will be concerned only with the small \x\ behaviour in contrast with the potential V\. Since the 
motion is semi-classical for small |x|, we expand the momentum p(x) as 




p(x) = J2m[E 2 + ^ « * *±+ Ehx (32) 



Notice the similarity between Eqn. ( pL 9[ ) and Eqn. (p2|). 

We will calculate the over the barrier reflection coefficient with E 2 > and small, which will be of interest later. 
Using the expansion in Eqn. (|32]), the semi-classical wave functions with the following boundary conditions, namely, 
that for x > the wave function is a right moving travelling wave while it is a superposition of an incident wave of 
unit amplitude and reflected wave for x < 0, are 

i>R = C 3 r 2+ ^« 2 / 2 (£>0) (33) 

ipL = (-zy^+h-* 2 / 2 + c 2 (~0 lE2+l2 e ie/2 (£ < 0) (34) 



where we have made the definitions 



mE$\ 1/4 V2mg 2 

-2g7) X ; e * = ^K~ (35) 



Following the variation of Eqn. (|3^) around an small semi-circle of radius p < \J g 2 j\E 2 \ (in contrast to the potential 
V\ where the radius p was large) in the upper half complex plane, we obtain, 

C 2 = C 3 exp^-7r£ 2 + ^ (36) 

Setting T — |C 3 | 2 = 1 and R = |C 2 | 2 , wc finally obtain, 

R = Te" 27rE2 = Te"^ 27rv/2 ™^ (37) 
Using the normalisation condition R + T = 1, we obtain, 

1 e -^2irV2mg 2 

T = , ; and R = = ^ (38) 

I _|_ e -^27rV2mg2 J _|_ e --jr 2 TT^ 2 mg 2 

Notice that the above result is independent of the energy E 2 and hence holds for E 2 < too. For small \x\, the lack 
of dependence on E 2 is not too surprising since the contour is such that it is not too close to the real turning points 
go = ± y / g 2 J\E 2 \. Anyway, when E 2 ~ + , p is "large" and therefore the contour is chosen to lie in the upper complex 
plane for the same reason as given in the analysis of the potential V\ in the previous section. 

Now, we will derive the result in Eqn. ( |37| ) using Eqn. The complex turning points qo are the roots of the 

equation E 2 = ~g 2 /q 2 ) where E 2 > and therefore, the turning points are go = ^\Zd2/E 2 = ±ip - Hence, we have 
to evaluate the integral, 
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J pdx = v / 2mE 2 J y 1 + ^dx (39) 

where the contour C encircles the point x — ipo in the upper half complex plane. However, since there is a singularity at 
the origin, we cannot deform the contour as was done when deriving Eqn. (|l4|). Therefore, as a means of regularisation, 
we modify the potential to 

V mod (x) = ^—r (40) 

x z + e z 

where the limit e — ► must be taken at the end of the calculation. The turning points for the modified potential are 
a^mod = ±i-\/e 2 + g 2 /E 2 while the poles of the modified potential are at x = ±ie < x mo d- Even in this case, there is 
a singularity on the path of integration which contributes to the integral rather than the turning point. Therefore, 
integrating upto +ie using the modified potential and back, we obtain, 



[ pdx = Hm2j2mE 2 [ \ 1+ f° dx 
Jc io f x + 6 



= Hm2iy/2mE 2 e dy^Jl + j ^ 
- 1 dy 



u 2 



2i^2mE 2 po 



— in \j l 2mE 2 pa = iir\j2mg 2 (41) 

We therefore recover the result given in Eqn. (|37|). From the above calculation it is clear that, due to the singularity 
at the origin, the reflection coefficient has no contribution from the turning point at all. 



III. PARTICLE PRODUCTION IN AN UNIFORM ELECTRIC FIELD 



We will now study a system consisting of a minimally coupled scalar field $ propagating in flat spacetime in an 
uniform electric field background. We consider two gauges, one a time dependent gauge while the other is a space 
dependent gauge and show how the tunnelling interpretation can be used to obtain the standard gauge invariant 
result in each case. We will not derive the standard result here. 



A. Time dependent gauge 

The four vector potential A 11 giving rise to a constant electric field in the x direction is assumed to be of the form 

^ = (0,-^,0,0) (42) 

The electric field is E = Eqx.. The minimally coupled scalar field $ propagating in flat spacetime, satisfies the 
Klein-Gordon equation, 

+ + iqA»)+m 2 )<S> = (43) 

where m is the mass and q is the charge of the field. The mode functions of <I> can be expressed in the form 
x) = /k(t)e x where fk(t) satisfies the equation, 

d 2 

^2 & + [™ 2 + k l + ( fc - + 1 E ^f] & = ! k ^ = (ky^z)- (44) 

Introducing the variables, 

t = ^q~Eht + (u/y/qEo) ; A = {k\ + m 2 )/qE Q (45) 

we obtain the equation, 
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-^/k-r 2 /k = A/ k (46) 

The above equation is essentially a Schrodinger equation in an inverted oscillator potential with a positive "energy" 
A. Since the energy is positive, the problem is essentially an over the barrier reflection problem. Using the results of 
section ( II A 2| ), we can calculate the reflection and transmission coefficients exactly as 



R = T~ x 5 T = --, " — r (47) 

I _|_ e -7 ™ 1 -J- 



where we have put % = 2m = ojq/2 = 1 and set a 2 = A in Eqns. ( |25| , 26). To identify the Bogoliubov coefficients a\ 
and (3\, we recast the normalisation condition R + T = 1 in the form, 

i-| = l (48) 
and then identify |/3a| 2 with R/T and |q;a| 2 with l/T. Therefore, the Bogoliubov coefficients are given by, 

|^|2 _ e -7rA _ e -Tr(k 2 L +m 2 )/qE 

\a X \ 2 = e-* X + 1 = e -Mk 2 ± +m 2 )/ q E„ + L ( 4Q ) 

Now, the transmission and reflection coefficients are time reversal invariant. They are dependent only on the energy 
(magnitude and sign) and for symmetric potentials, independent of the direction in which the boundary conditions 
are applied. To obtain a dynamical picture of particle production, we have to interpret these quantities suitably. 
In the present case, the following interpretation seems adequate. A purely positive frequency wave with amplitude 
square T in the infinite past, t — > — oo, evolves into a combination of positive and negative frequency waves in the 
infinite future t — > oo with the negative frequency waves having an amplitude square R and the positive frequency 
waves having an amplitude unity. The quantity R/T determines the overlap between the negative frequency modes in 
the distant future and the positive frequency modes in the distant past (the notation here differs from the treatment 
given in || , . This is identified with the modulus square of the Bogoliubov coefficient j3\ which is the particle 
production per mode A. Using the normalisation condition satisfied by the Bogoliubov coefficients, \a\\ 2 — |/3a| 2 = 1, 
|q:a | 2 can be calculated to be l/T. Once the Bogoliubov coefficients have been identified, the effective lagrangian can 
be easily calculated. This derivation will not be repeated here. We refer the reader to Ref. and Ref. Q for the 
explicit calculation. 

Note that the particular interpretation given in this case is due to its similarity with the more rigourous calculation 
by quantum field theoretic methods. In the next section, where we discuss the space dependent gauge, we will be 
forced to adopt a different interpretation in order to identify particle production. 

B. Space dependent gauge 

The four vector potential giving rise to a constant electric field in the x direction is now assumed to be of the 
form 

A* = (-E x, 0,0,0) (50) 
The electric field is E = EqX.. The field $ satisfies Eqn. ( pf3| ) as before. Substituting for the potential v4 M from 



Eqn. (50) into Eqn. (|4^), we obtain, 

(<9 t 2 - V 2 - 2iqE xd t - q 2 E$x 2 + to 2 ) $ = (51) 

We write $ in the form 

$ = e- lu)t e ik * y+lk * z <t>{x) (52) 
and obtain the differential equation satisfied by (j) as 

^ + ((w + qE x) 2 -k 2 ± - m 2 ) = (53) 
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where we have used the notation k\ = k 2 + k 2 . Making the following change of variables 



p = x/^x + {u/y/qE^j ; A — (kj_ + m 2 )/qE (54) 
into the differential equation for 4>, it reduces to the form, 

~%- p ^ = ~ M (55) 
In this form, we see that the above differential equation has the form of an effective Schrodinger equat ion with an 



inverted harmonic oscillator potential and an negative energy — A. If we apply the results of section ( II A 2 ), we obtain 
the result for tunneling through the barrier. Following the treatment in Ref. Q and using the results of section ( II A 2| ) 
we can calculate the reflection and transmission coefficients exactly as 



R = r ; T = r (56) 

where we have put fi — 2m = ljq/2 = 1 and set a 2 — A in Eqns. ( |25| , |26| ). We cast the renormalisation condition 
R + T = 1 in the form 

1 T 

R~ R =1 < 57 > 

and then identify the rate of particle production per mode with T/R. The interpretation of particle production using 
the tunnelling picture now proceeds as follows. A right moving travelling wave of amplitude square 1/ R is incident on 
the potential. A fraction T/R is transmitted through it and a wave of unit amplitude is scattered back. The tunnelling 
probability, which is T/R, is interpreted as the rate at which particles are being produced by the background electric 
field. This matches exactly with the expression for \/3\\ 2 given in Eqn. (|49|). With this interpretation, we recover the 
usual gauge independent result. 

In summary, it is seen that by a judicious choice of interpretation of the transmission and reflection coefficients in 
each of the two gauges, the standard gauge invariant result can be obtained. In fact, the tunnelling interpretation 
can also be used for gauges of the form 

A p = {—Eqx ± E\t, 0, 0, 0) ; = (0, —E t ± Eix, 0, 0) (58) 

where the condition |Ei| ^ \Eq\ holds strictly. Depending on the magnitude of E\, the problem reduces to either an 
"over the barrier reflection" or a "tunnelling through the barrier" process. The transmission and reflection coefficients 
do not depend on E\ as it ought to be. But the tunnelling interpretation is not always valid. For instance, setting 
| .Ei | = j-Eol in the above gauges reduces the problem to a first order differential equation instead of a second order one. 
(This case will be considered later in a future publication.) Further, one can show that even in a time independent 
magnetic field, solving the effective Schrodinger equation gives non-zero transmission and reflection coefficients 
Since these coefficients are non-zero, applying the tunnelling picture implies particle production which is contrary to 
the result obtained by Schwinger in jl|. Hence, the tunnelling picture does not always produce consistent results. 
In the next section, we look at Schwarzchild-like spacetimes with a horizon. Here too, we construct the tunnelling 
picture near the horizon in order to explain particle production. We shall see that such a picture does not produce 
the correct results. 



IV. PARTICLE PRODUCTION IN SPACETIMES WITH HORIZON 



Hawking's result that a black hole radiates is essentially a semi-classical result. The thermal radiation results 
because of the presence of a horizon in the spacetime structure. We will review briefly the conventional derivation 
of the thermal radiation using path integrals. Consider a patch of spacetime, which in a suitable co-ordinate system, 
has one of the following forms (we assume c = 1): 

ds 2 = B{r)dt 2 - B-\r)dr 2 - r 2 (d8 2 + sm 2 {6)d<j) 2 ) (59) 

or 

ds 2 = B{x)dt 2 - B-\x)dx 2 - dy 2 - dz 2 (60) 
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where B(r) and B(x) are functions of r and x respectively. The horizon in the above spacetimes is indicated by the 
surface r — r$ (x = xq) at which B(r) (B(x)) vanishes. We further assume that B'(r) — dB/dr (B'(x) = dB/dx) is 
finite and non-zero at the horizon. Co-ordinate systems of the form ( |59| ) can be introduced in parts of the Schwarzchild 
and de-Sitter spacetimes while that of the form (|60|) with the choice B{x) = 1 + 2gx represents a Rindler frame in 
flat spacetime. Given the co-ordinate system of ( p9[ ) say, in some region R, we first verify that there is no physical 
singularity at the horizon, which in the case of the Schwarzchild black hole, is at the co-ordinate value tq = 2M where 
M is the mass of the black hole. Having done that, we extend the geodesies into the past and future and arrive at 



two further regions of the manifold not originally covered by the co-ordinate system in (59) (the Kruskal extension). 
It is now possible to show that the probability for a particle with energy E to be lost from the region R in relation 
to the probability for a particle with energy E to be gained by the region R is given by the relation 

P. - P e~$ E fBll 

J loss — J gain c V u / 

where (3 = 8irM. This is equivalent to assuming that the region R is bathed in radiation at temperature In the 
derivation given in the paper by Hartle and Hawking Q , thermal radiation is derived using the semiclassical kernel by 
an analytic continuation in the time co-ordinate t to complex values and it was shown that the probability of emission 
(loss) from the past horizon was related to absorption (gain) into the future horizon by the relation ([n]) . 

Since all the physics is contained in the (t, r) plane (or the (t, x) plane), we will discuss Hawking radiation in 1 + 1 
dimensions in the following sections and later show in the appendix that the results generalise naturally to 3 + 1 
dimensions. We first derive the semi-classical result in the (r, t) (or (x,t)) plane by applying a certain prescription 
to bypass the singularity encountered at the horizon. After this, we reduce the problem of the Klein-Gordon field 
propagating in the Schwarzchild spacetime to an effective Schrodinger problem in (1+1) dimensions and rederive 
the semi-classical result by using the same prescription on the semi-classical wave functions. Then the results of 
section (jllj) are used to compute the ratio of the transmission and reflection coefficients and we show that this ratio 
satisfies an equation of the form 

R = Te-^ E (62) 

where /3 _1 = (87rM) _1 is the expected Hawking temperature. This can be interpreted to mean thermal radiation in 
exactly the same way as Eqn. ( |6l| ) but with a temperature twice the expected value. This proves that the tunnelling 
picture, when naively applied, does not give the same result as the semi-classical analysis does. A possible qualitative 
explanation of this result will also be given. 



A. Hawking Radiation in 1+1 dimensions 

Consider a certain patch of spacetime in (1+1) dimensions which in a suitable co-ordinate system has the line 
element, (with c = 1) 

ds 2 = +B(r)dt 2 - B- l (r)dr 2 (63) 

where B(r) is an arbitrary function of r. We assume that the function B(r) vanishes at some r = ro with B'(r) 
fill fir being finite and nonzero at r$. The point r — tq indicates the presence of a horizon. It can be easily 
verified that no physical singularity exists at the horizon since the only component of the curvature tensor IZtrtr = 
— (l/2)(d 2 B(r)/dr 2 ) does not become infinite at the horizon. Therefore, near the horizon, we may expand B(r) as 

B(r) = B'(r )(r - r ) + 0[(r - r ) 2 } = R(r )(r - r ). (64) 

We now consider a minimally coupled scalar field $ with mass too propagating in the spacetime represented by the 
metric (^33|). The equation satisfied by the scalar field is, 

□ + ^)$ = (65) 

where the □ operator is to be evaluated using metric (|63|). Expanding the LHS of equation J65|) , one obtains, 

1 <9 2 $ d ( , N <9$\ m 2 , , 

l B (r)ir- =-TT $ - (66) 



B(r) dt 2 Or \ Or 
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The semiclassical wave functions satisfying the above equation are obtained by making the standard ansatz for $ 
which is, 



$(r, t) = exp 



- h S(r,t) 



(67) 



where S is a function which will be expanded in powers of h. Substituting into the wave equation (p6|), we obtain, 



1 



dS 



B(r) \ dt 



B(r) 



dS_ 

dr 



1 d 2 S _ cPS_ _ dB{r) dS 



= 



Expanding S in a power series of (h/i), 



S(r,t) = S {r,t) + - Sx(r,t) + 



S 2 (r,t). 



and substituting into Eqn (^8|) and neglecting terms of order (h/i) and greater, we find to the lowest order, 

2 / no N 2 



1 



dSo 
B(r) \ dt 



-B(r) 



(dSoY a 



(68) 



(69) 



(70) 



Eqn (70) is just the Hamilton- Jacobi equation satisfied by a particle of mass too moving in the spacetime determined 
by the metric (|63|). The solution to the above equation is 



S (r,t) = -Et± ^ ' ^WfP-mgBtr) 



(71) 



where E is a constant and is identified with the energy. Notice that in the case of too = 0, Eqn (70) can be exactly 
solved with the solution 



S mo =o(r : t) =F 1 (t- r*) + F 2 (t + r*) 
where the "tortoise" co-ordinate r* is defined by 

dr 



B( r y 



(72) 



(73) 



and Fi and F 2 are arbitrary functions. If F\ is chosen to be F\ = —Ft + Er* and F 2 chosen to be F 2 = —Ft — Er*, 
then it is clear that the solution given in Eqn (72) is the same as that in Eqn ( |7l| ) with too set to zero. Therefore, in 
the case too = 0, the semiclassical ansatz is exact. In the following analysis we will specialise to the case mo = for 
simplicity. The case too 7^ will be considered later. The essential results do not change in any way. 

Knowing the semiclassical wave function, the semiclassical kernel K(r 2 , t 2 \ 7*1, t\) for the particle to propagate from 
(ti, r-y) to (t2, r 2 ) in the saddle point approximation can be written down immediately. 



(74) 



K(r 2 ,t 2 \r x M) = TV exp ( -S^M'^xM , 



where S is the action functional satisfying the classical Hamilton- Jacobi equation in the massless limit and TV is a 
suitable normalization constant. S(r 2 ,t 2 ;ri,ti) is given by the relation 



S(r 2 ,t 2 ;ri,h) = S(2, 1) = -E(t 2 - h) ± E 

Jn 



dr 

WY 



(75) 



The sign ambiguity (of the square root) is related to the "outgoing" (dS/dr > 0) or "ingoing" (dS/dr < 0) nature 
of the particle. As long as points 1 and 2, between which the transition amplitude is calculated, are on the same side 
of the horizon (i.e. both are in the region r > ro or in the region r < tq), the integral in the action is well defined 
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and real. But if the points are located on opposite sides of the horizon then the integral does not exist due to the 
divergence of P _1 (r) at r = ro- 

Therefore, in order to obtain the probability amplitude for crossing the horizon we have to give an extra prescription 
for evaluating the integral pj. Since the point B = is null we may carry out the calculation in Euclidean space or 

cquivalcntly — use an appropriate ie prescription to specify the complex contour over which the integral has to be 
performed around r = ro . The prescription we use is that we should take the contour for defining the integral to be an 
infinitesimal semicircle above the pole at r = ro for outgoing particles on the left of the horizon and ingoing particles 
on the right. Similarly, for ingoing particles on the left and outgoing particles on the right of the horizon (which 
corresponds to a time reversed situation of the previous cases) the contour should be an an infinitismal semicircle 
below the pole at r — ro- Equivalently, this amounts to pushing the singularity at r — ro to r — ro =p ie where the 
upper sign should be chosen for outgoing particles on the left and ingoing particles on the right while the lower sign 
should be chosen for ingoing particles on the left and outgoing particles on the right. For the Schwarzchild case, this 
amounts to adding an imaginary part to the mass since ro = 2M. 

Consider therefore, an outgoing particle (dS/dr > 0) at r = r\ < ro- The modulus square of the amplitude for 
this particle to cross the horizon gives the probability of emission of the particle. The contribution to S in the ranges 
( r i 5 r o ~ e) and (ro + e, r-i) is real. Therefore, choosing the contour to lie in the upper complex plane, 



f r ° +e dr 

Remission] = — lim / „. . + (real part) 
^oj ro _ e B(r) 



ittE 

= —. — r + (real part) (76) 
R[r ) 

where the minus sign in front of the integral corresponds to the initial condition that dS/dr > at r = r± < ro. For 
the sake of definiteness we have assumed P(ro) in Eqn (|64| ) to be positive, so that B(r) < when r < ro. (For the 
case when R < 0, the answer has to be modified by a sign change.) The same result is obtained when an ingoing 
particle (dS/dr < 0) is considered at r = r\ < ro- The contour for this case must be chosen to lie in the lower 
complex plane. The amplitude for this particle to cross the horizon is the same as that of the outgoing particle due 
to the time reversal invariance symmetry obeyed by the system. 

Consider next, an ingoing particle (dS/dr < 0) at r = r 2 > r . The modulus square of the amplitude for this 
particle to cross the horizon gives the probability of absorption of the particle into the horizon. Choosing the contour 
to lie in the upper complex plane, we get, 



[ r °~ e dr 

S [absorption] = — lim / + (real part) 

e^oj +e B{r) 



= -- — - + real part 77 
R[ro) 

The same result is obtained when an outgoing particle (dS/dr > 0) is considered at r — r^ > ro- The contour for 
this case should be in the lower complex plane and the amplitude for this particle to cross the horizon is the same as 
that of the ingoing particle due to time reversal invariance. 

Taking the modulus square to obtain the probability P, we get, 

f 2tt£ \ , x 

P [emission] cx exp —-^7 — r ] (78) 



and 



so that 



R{ro) 



P [absorption] cx exp ( ) (79) 



R(ro) 



( AttE \ 

P [emission] = exp — - P [absorption]. (80) 

V R( r o)J 

Now time reversal invariance implies that the probability for the emission process is the same as that for the absorption 
process proceeding backwards in time and vice versa. Therefore we must interpret the above result as saying that the 
probability of emission of particles is not the same as the probability of absorption of particles. In other words, if the 
horizon emits particles at some time with a certain emission probability, the probability of absorption of particles at 
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the same time is different from the emission probability. This result shows that it is more likely for a particular region 
to gain particles than lose them. Further, the exponential dependence on the energy allows one to give a 'thermal' 
interpretation to this result. In a system with a temperature the absorption and emission probabilities are related 
by 

P[emission] = exp(—/3E)P [absorption] (81) 

Comparing Eqn (^TJ) and Eqn (|8(i|), we identify the temperature of the horizon in terms of R(ro). Eqn ( |80| ) is based 
on the assumption that R > 0. If R < there will be a change of sign in the equation. Incorporating both the cases, 
the general formula for the horizon temperature is 

r * = § (82) 

For the Schwarzchild black hole, 

B{r) = (l - « JL(r 2M) + 0[(r - 2M) 2 ] (83) 

giving R = (2M) _1 , and the temperature /3 -1 = 1/8ttM. For the de-Sitter spacetime, 

B(r) = (1 - H 2 r 2 ) w 2iJ(iJ" 1 - r) = -2#(r - If -1 ) (84) 

giving P^ 1 — H/2tt. Similarly for the Rindler spacetime 

B(r) = (l + 2gr) = 2g(r+(2g)- 1 ) (85) 

giving = g/2ir. The formula for the temperature can be used for more complicated metrics as well and gives the 
same results as obtained by more detailed methods. 

The prescription given for handling the singularity is analogous to the analytic continuation in time proposed by 
Hawking Q| to derive Black hole radiance. If one started out on the left of the horizon and went around the singularity 
r = ro by a 2tt rotation instead of a rotation by n, it can be easily shown that it has the effect of taking the Kruskal 
co-ordinates (v,u) to (—v,—u). A full rotation by 2ir around the singularity can be split up into two parts to give 
the amplitude for emission and subsequent absorption of a particle with energy E. Since the amplitudes for the two 
processes are not the same in the presence of a horizon, one obtains the usual Hawking radiation given in Eqn (p0|) 
with the value of R(r ) being (2M) _1 . This process is similar to that given in |4f| which relates the amplitudes 
involving the past and future horizons. In Hawking's paper, analytically continuing the time variable t to t — AMin 
takes the Kruskal co-ordinates (v, u) to (— v, — u) and since the path integral kernel is analytic in a strip of AMiir 
below the real t axis, Hawking radiation is obtained by deforming the contour of integration appropriately. 

When too =^ 0, the validity of the semi-classical ansatz must be verified. To do this, consider the perturbative 
expansion (|69|). Retaining the terms of order Ti/i and neglecting higher order terms, one finds, upon substituting for 
So given by the relation (|7l]) and solving for Si, 

S, = -*« ± EEl J ^r ) 7W ^ m - \ HE' - mSB(r)) (86) 

where E\ is a constant. From the above equation, it is seen that Si has a singularity of the same order as So at r = tq. 
When calculating the amplitude to cross the horizon, the contribution from the singular term just appears as a phase 
factor multiplying the semiclassical kernel and is inconsequential. The non-singular finite terms do contribute to the 
kernel but they contribute the same amount to S [emission] and S [absorption] and they do not affect the relation 
between the probabilities P[emission] and P [absorption]. Subsequent calculation of the terms S2, S3, and so on, show 
that all these terms have a singularity at the horizon of the same order as that of So- Their contribution to the 
probability amplitude is just a set of terms multiplied by powers of h which can be neglected. From this we may 
conclude that the semiclassical ansatz, in the perturbative limit, is a valid one. 



B. Reduction to an effective Schrodinger Problem in (1 + 1) dimensions 

Consider the relativistic equation for the wave function <!> in Eqn ( |66| ) . We include the mass Too here but later we 
shall see that it does not appear in the final reduced Schrodinger equation. We first set $ = e~ lEt ^ h ^{r) to obtain, 
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B(r) 



d 2 ^ dB(r) 



dr 2 



dr dr 



E 2 



ml 



We then make the substitution 



fTB(r) h 
*(r) = P(r)Q(r) with P(r 



* = 



(87) 



to get the equation, 



d 2 Q{r) 
dr 2 



B"(r) (B'(r)f 



E 2 



2B(r) AB 2 {r) h 2 B 2 (r) h 2 B(r) 



Q(r) = 



(89) 



where B' denotes dB /dr and B" denotes d 2 B/dr 2 . Near the horizon, we use the expansion of B(r) given in Eqn (]64|). 
Neglecting terms of order l/(r — r ) as compared to terms of order l/(r — r ) 2 , we get, 



d 2 Q(r) 
dr 2 



(r - r ) ; 



rQ(r) = 



where 



.9 



£ 2 



^ 2 i? 2 



(90) 



Notice that mo appears nowhere in the above equation. Very close to the horizon, the term containing the mass does 
not contribute significantly. Eqn (|90| ) is therefore applicable to both massless and massive scalar particles. Making 
the simple transformation x — (r — r ), we finally obtain the effective Schrodinger equation for the system with 
H = 2m = 1 with a potential {—g/x 2 ). 



d 2 Q(x) g 
dx 2 x 2 



Q{x) = 



(91) 



This potential is symmetric about the origin x — 0. We now include an "energy" E (not to be confused with the 
energy E of the field in the original relativistic system) to obtain, 



d 2 Q(x) 
dx 2 
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Q(x) = EQ(x). 



(92) 



The energy E, which can be positive or negative, has been included in order to consider fully the properties of the 
potential —g/x 2 . We will finally be interested in the case E = which is the case of interest. (The energy spectrum 
is continuous for all values of E which, for E < 0, is surprising since for energies less than the potential energy, the 
spectrum is usually discrete.) In the subsections to follow we will rederive the semiclassical result giv en in th e previous 
section and then calculate the reflection and transmission coefficients using the results of section (II A 3). We then 
look for a correspondence between the reflection and transmission coefficients with the semiclassical results. 



1. Semiclassical analysis of the effective Schrodinger equation 



The semiclassical analysis follows closely the method adopted in section ( I V A ) . The action functional A for a 
classical particle moving in a potential —g/x 2 satisfies the Hamilton- Jacobi equation 



dA fdA\ 2 g 

m + 

The solution can be immediately written down as 



dt ' \ dx J x 2 ^ 



A = -Et± J ^-^Ex 2 +g (94) 

Eqn (^) has an integral which is divergent if the action is computed for points lying on the opposite sides of the horizon 
x = 0. Since this has a similar form to Eqn (|7l"|), the prescription used in evaluating ,!? [emission] and S [absorption] 
can be similarly used to evaluate ^[emission] and ^[absorption]. The results are 

^[emission] = iTty/g + (real part) 
^[absorption] = —in^/g+ (real part). (95) 
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Constructing the semiclassical propagator as before and taking the modulus square to obtain the probabilities for 
outgoing and ingoing particles, we get 



P [emission] = exp 



-47T 



E 2 



fi 2 R 2 



P [absorption] 



(96) 



Note that the result is independent of E. Therefore, the above result holds even in the case E = 0. The formula 
above does not represent a thermal system since the energy E appears in the squareroot. For large energies such that 
(E 2 /h 2 R 2 ) » (1/4), the above equation can be written in the form that suggests a thermal system 



P [outgoing] = exp 



A-kE 



fiR 



P [ingoing] 



(97) 



The temperature for the system is the same as that in Eqn fl82[ ) and one recovers the usual result. For the case 
of the Schwarzchild black hole, the condition [E 2 /h 2 R 2 ) (1/4) reduces to {E/h) ^ l/AM. In proper units, it can 
be written as 



he 3 
AGM 



= M P & 



,Mp 
AM 



(98) 



where Mp — yJfic/G is the Planck mass. Since the energy E must be much smaller than the Planck energy Mpc 2 in 
the semiclassical limit, the condition above reduces to M 3> Mp i.e. the mass of the black hole must be far greater 
than the Planck mass for the spectrum to be thermal. This differs significantly from the result given in Eqn ( |80| ) 
which does not put any condition on the size of the black hole for the spectrum to be thermal. To verify that the 
semi-classical analysis is valid, one must compute the correction terms and check that these have a singularity of the 
same order as possessed by A. To do this, consider the effective Schrodinger equation (p2|) with factors of h put in. 



_ 2 d 2 Q{x) 



4<2(z) = EQ(x) 



dx 2 x 

Putting Q(x) = exp(iA(x)/h), and substituting into Eqn (B9), 



-%h 



d 2 A(x) 
dx 2 



dA(x 
dx 



= E 



Expanding A in powers oifi/i, we get, 



A = A+ -A x 



A 2 + 



(99) 



(100) 



(101) 



Substituting into Eqn (100) and proceeding as usual, we find that A is given by Eqn (94). The next term Ai is given 
by 



Ax =g 



dx 



1 



x Ex 2 + g 



(102) 



The relation for Ax also has a singularity at the origin of the same order as A. Explicit calculation of the subsequent 
terms in the expansion of A reveals that all these terms have a singularity of the same order as that of A and therefore 
their net contribution to the kernel is either as phase factors or as the exponential of finite terms multiplied by powers 
of h. Therefore, we conclude as before that the semiclassical ansatz is valid. 



2. Tunneling Interpretation 



Instead of the purely semiclassical derivation given above, we attempt to see if the relation between the reflection and 
transmission coefficients for the eff ective potential calculated above could be interpreted to give a thermal spectrum. 
Using the results in section ([I A3) and in particular Eqn. (p7|) after setting 2m = h = 1, we obtain, 



R = Te- 27T ^ = Texp 



2ttE 

Yr 



(103) 
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for "large" energies satisfying the inequality Eqn. ( |9q ) . We may interpret the transmission coefficient as the probability 
of emission of particles while the reflection coefficient as that for absorption but we obtain a temperature which is 
twice the temperature obtained using the semiclassical analysis. It is also clear why this is so. The transmission and 
absorption of particles to and from the horizon is clearly not symmetric. In the semiclassical analysis, this asymmetry 
was put in by hand in the form of the prescription used to evaluate the semiclassical kernel while in the tunneling 
analysis, the potential near the horizon is completely symmetric and there is no apriori reason for the transmission 
coefficient to be different on cither side of the horizon. The tunnelling interpretation, in the case of the electric field, 
works in both gauges because of the structural similarity of the effective Schrodinger equation in both cases. Further, 
in the time dependent gauge, boundary conditions play an important role in determining particle creation. The 
boundary condition that, at t — — oo, the quantum field is initially in a vacuum state, introduces an asymmetry into 
the problem which gives rise to particle production in the infinite future as t — > oo. But, in the space dependent gauge, 
such a time asymmetry in the boundary conditions cannot be implemented since the problem is time independent. 
But, since the structure of the effective schrodinger equation to be solved is the same as that in the time dependent 
gauge, suitable ratios of the transmission and reflection coefficients, which are symmetric in both space directions, 
can be identified with the Bogolibov coefficients with the subsequent recovery of particle production. 

In the Schwarzchild spacetime, the analogue of the "time dependent gauge" is provided by the study of a collapsing 
shell of matter [0 which, close to the black hole horizon, gives rise to production of particles which have a thermal 
distribution in energy. In this case, the boundary condition considered at the collapsing shell causes an exponential 
redshift of incoming modes. The Bogolibov coefficients can be calculated in a straight forward manner and \f3\\ 2 is 
Planckian in form and particle production is established in the same manner as was done in the electric field case. The 
analogue of the "space-dependent gauge" is provided by the system where an eternal black hole is in equilibrium with 
a bath of thermal radiation. The transmission and reflection coefficients are easily determined in this case. B ut, i f 
we attempt to identify the Bogoliubov coefficients using ratios of these coefficients, we obtain the result in Eqn. ( |103| ) 
where the tun nelling temperature is twice Hawking's result. On the other hand, if we use the correspondence given 
in section (III) for the identification we obtain results similar to those obtained for the electric field and this does not 
yield a thermal spectrum for the emitted particles. A possible qualitative explanation for the factor 2 in the tunnelling 
tempe ratu re is that, since there are two disjoint Schwarzchild sectors in the full Kruskal manifold, the energy E in 
Eqn. (103) should be multiplied by a factor 2 resulting in the recovery of Hawking's result. Therefore, the radiation 
energy seen at infinity is the result of tunnelling in both sectors. Such an explanation holds for all the spacetimes 
considered in this paper with a horizon, namely, the Rindler and the de-Sitter where the full manifold consists of 
two mutually disjoint but identical spacetimes. But a more quantitative and satisfactory tunnelling interpretation in 
black hole like spacetimes is needed. 



V. COMPLEX PATHS INTERPRETATION OF PARTICLE PRODUCTION IN ELECTRIC FIELD 



In section (III), we noted how particle production can be calculated using the tunnelling interpretation. This 
interpretation gives the same result in both the space as well as the time dependent gauges. The spectrum of particles 
produced by the electric field is not thermal in contrast to the spectrum seen by a Rindler observer. We use the 
formal tunnelling method to show, in a heuristic manner, how this particle production can be obtained by tunnelling 
between the two sectors of the Rindler spacetime. 

The imaginary part of the effective lagrangian ImL c ff , which is related to the probability of the system to remain 
in the vacuum state for all time is given by |^] 

v 2 * 1 (^o) 2 (-1)" +1 ( 7rm 2 \ , , 

ImL eff = J2o to % 2 ex P 104 

where m is the mass, q is the charge and E is the magni tude of the electric field. We will derive the above expression 



for ImLcfj using the general arguments given in section (II A) 



Consider the Hamilton- Jacobi equation for the motion of a particle in an electromagnetic field in (1 + 1) dimensions 
with a proper time co-ordinate s. 

1 fdF At \ 2 1 fdF \ 2 1 2 dF . . 

2-U+^ J ~2 -2™ +^=° (1 ° 5) 

where F is the action and = (A 1 , A x , 0, 0) is the four vector potential. We have neglected the dependence of F on 
the y and z co-ordinates. This will be justified later on. In a time-dependent gauge, given by Eqn. (fl2|), the action F 
can be easily solved for by using the ansatz F = —Es + p x x + f(t) to give, 
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ftl 

F{t 1) x 1 ;t ,x ;s)=-Es+p x {x 1 -x )± dt y/{p x + qE t) 2 + (m 2 + 2E) (106) 

J to 

where p x is the momentum of the particle in the x direction and E is the energy of the particle corresponding to the 
proper time co-ordinate s. The trajectory of the particle in the (t, x) plane is the usual hyperbolic trajectory given 
by, 

(t - U) 2 -{x- x,f = -fcg (107) 

where ti, Xi and ko are suitable constants. For any fixed position x, there are thus two disjoint trajectories corre- 
sponding to motion in the two Rindler wedges. 

Let us consider the tunnelling of a particle from one Rindler trajectory to the other and back in the imaginary time 
co-ordinate t. This means that the particle comes back to the same spacetime point as it started from. Sin ce the 
process is a tunnelling process, the proper time taken is zero. Therefore, choosing the positive sign in Eqn. (|l06| ) (this 
choice give a tunnelling probability that is exponentially damped), we have, 

m 2 + 2E 



F{t Q ,x o ;t o ,x o ;0) = j dt^(p x + qE t) 2 + (m 2 + 2E) = — j> du^l + u 2 
i q> dr\/l - t 2 = i / ddcos 2 (9) 



qE J qE 
in(m 2 + 2E) 
qEo 



(108) 



where we have made the following changes of variable (p w + qE^t) = u = it and r = sin(0). Taking the limit 
E — > 0, the expression for exjp(iF) from the above equation is seen to be exactly the same as the exponential term in 
Eqn. (104) for n = 1. The same argument can be repeated for the particle tunnelling n times to and fro to give 

F n (t ,x o ;t ,x Q ;0) = i 1 ™^E d9cos 2 {9) = m{m2 + 2E) n. (109) 
qE Jo qE Q 



Again, in the limit E — > 0, exp(iF n ) is seen to match with the exponential part of the nth term in Eqn. (104). 
Therefore, the imaginary part of the total effective lagrangian can be written down immediately as 



ImLeff = (prefactor) exp ( — n ) (HO) 

^ v qE ° / 

where the prefactor can only be c alculated using the exact kernel. The dependence of the prefactor on n and the 
phase factor (— l) ,l+1 in Eqn. (104) can be deduced using the following arguments. 

The formal expression of the path integral kernel for the above electric field problem in the time dependent gauge 
is given by |g 

K(a,b;s) = (a\e ish \b) (111) 

where K(a,b;s) is the kernel for the particle to propagate between the spacetime points a = {x , x) and b = (y , y) 
in a proper time s and h is the hamiltonian given by 

h = ^(idi - qAiXid* - qA*) - ^m 2 (112) 



where A 1 is the four vector potential given in Eqn. (|42j) and q and m are the charge and mass of the particle respectively. 
Going over to momentum co-ordinates and considering the coincidence limit x = y, the kernel can be written in the 
form, 



K{x\y Q -^s) = --^ r^ gi x ,y ;s) (113) 



where Q(x°,y°; s) is given by 

g(x ,y ;s) = {x°\e isH \y ) (114) 
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1 I u ' „.x9 2 \ / ,2 7^2 2, 2 



and H is the hamiltonian 

H = ~ 2 + iPx + qEat)1 + "'"J = " 2 [dp 2 " + qAE ° pZ + (115) 

with p — (t + p x /qEQ). In the expression for the kernel in momentum co-ordinates, we have integrated over the 
transverse momentum variables p y and p z . The above hamiltonian is that of an inverted harmonic oscillator. Since 
all reference to p y and p z have disappeared in H , the dependence of F on t he y and z co-ordinates was neglected 
when writing down the expression for the Hamilton- Jacobi equation in Eqn. ( |105| ). The expression for the effective 
lagrangian is then given by 



L c s = —i r-K(x°,x°; X ,x; S ) = --^ f°° ^ C dp x g(x°,x°; 
Jo s 2(2tt) 2 J s 2 



s) (116) 



We would like to evaluate the propagator Q(x°,x°; s) for a tunnelling situation where the particle tunnels from the 
point x° and back in loops. Since the path integral is not well defined for closed loops, it will have to be evaluated in 
some approximate limiting procedure which is outlined below. 



Since the tunnelling potential is that of an inverted oscillator, we can use all the results of section (|IA_2j). The 
semi-classical wave functions are given in Eqns. ( p0| , ^l| ). We would like to first account for the factor (—1)" that 
arises when a particle tunnels from one side of the barrier to the other and back. Consider an incident wave to the 
right of the barrier and impinging on it. Using the method of complex paths, we rotate this wave in the lower complex 
plane (this is the only route possible for the same reason as that given when rotating a right moving travelling wave 
in the upper complex plane) to obtain a wave again incident on the barrier with a energy independent phase factor 
exp(i7r/2) being picked up (other factors dependent on the energy are also picked up but these are not important 
here). Since this wave is moving in the wrong direction, we assume that the particle that has tunnelled through has 
the same amplitude as the rotated wave but is moving away from the barrier. This just involves changing the sign of 
the argument of the exponential in the expression for the rotated wave. Rotating this left moving wave again in the 
upper complex plane now, the final wave obtained is a right moving wave with another extra phase factor of exp(j7r/2) 
being picked up. The total phase change with respect to the incident wave is thus exp(«7r). Since this phase factor is 
independent of the energy, the propagator for the tunnelling process too, after one such rotation, will pick up a phase 
of exp(i7r). Similarly for n rotations, n taking the values 1, 2, 3, . . ., the phase acquired will be exp(in7r) = (—1)". 

Therefore, the propagator Q for n loops, Gn(x°, x°;s), can thus be written as 

G n (x°, x°; s) = N(p x , to, E )e™ e lF ^ x ° > x ° ^ (117) 

where F n (x° ,x°;s) is the classical action for the tunnelling process and N is the prefactor that arises after evaluating 
the "sum over paths". This prefactor is not expected to depend on the proper time s since the tunnelling process 
takes place instantaneously or on the number of rotations n. So the only q uan tities it may depend on are p x , m and 
E. Now the expression for the classical action for n loops is given in Eqn. (|109|) with E set to zero. Substitutin g thi s 



into the expression for the propagator in Eqn. (117) and evaluating L c g(n) for n loops from the relation in Eqn. (116) 
we obtain infinity as the answer which is meaningless. Therefore, we adopt the following limiting procedure in order 
to obtain a finite and meaningful answer. Consider the action for the tunnelling process. 

F = -Es +p x x± J dty/(p x + qE Q t) 2 + (to 2 + 2E) (118) 
Choosing the positive sign and setting p x + qE^t = iy/m 2 + 2E sin(#), we obtain, 

F = -Es + p x x ± i / d9(l + cos(20)) (119) 



2qE 

For closed paths only, with 9 taking the values from to 2n7r, the above action can be written as, 

F n = -Es±i r ^-^6 (120) 

We have thrown away the p x x term while retaining the —Es term since the dependence on the x co-ordinate is really 
irrelevant. Defining a new variable 6 = i9/2qEo and rescaling s = as' , the above action becomes 

F n = -Eas' + (to 2 + 2E)6. (121) 
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Choosing a appropriately, the above action can be cast into a form that matches the action for a fictitious free particle 
in (1 + 1) dimensions with "energy" aE and "momentum" (m 2 +2E) satisfying the classical energy-momentum relation, 

aE = i(m 2 + 2E) 2 (122) 

where the particle's "mass" is set to unity for convenience. The above equation determines the quantity a. Therefore, 
F n can also be written in the form 

Fn(M - i; , ) = <^ = «£ (123) 

where 9\ and (9 2 are the initial and final states of the free particle with s' being the p rope r time taken. (Note that 
(82 — Oi) = 2inir/2qEo.) Substituting this into the expression for Gn (x°, x°; s) in Eqn. (117) and evaluating only the 
integral over s in the expression for the effective lagrangian in Eqn. (|116|) without taking the limits, we obtain 



/ 



dS g n ( X °, A s) = N{pxj m> E yn, f ds ^ (ia{h - hf 



-N(p x , m, E)j™ , - 2 exp (iF n (6 2 , g i; a)) . (124) 
ia(0 2 ~ #1) 



Notice that the prefac tor t o the exponential has no dependence on the proper time s. Now, we use the form for the 
action given by Eqn. (120) and taking the limits for s from to 00, the effective lagrangian for n loops, L Q g(n) can 
be written in the form, 

, , * (qE ) 2 (-l) n+1 ( irm 2 \ 8 / 2ttE \ [°° d Px , x . . 

L eff(n) = ^-TTT^r- ^exp —n -exp -n / -^N{ Px ,m,E) (125) 



2 (2tt) 3 n 2 * V qE J a ^\ qE J J_ x 2ir 



where we have set (62 — Oi) = 2inn/2qEo. Taking the limit of E — >• and using the expression for a in Eqn. ( 122 ), 
we find that N must satisfy the relation 



lim - — s -77 f -!—-N(p T ,m,E) = 1 

B^o (m 2 + 2E) 2 2n Kl x ' ' ' 



(126) 



so that the imaginary part of the effective lagrangian for n loops, L e g(n), matches the nth term in Eqn. (104). 
Therefore, in this manner, the contributions to the imaginary part of the effective lagrangian for the uniform electric 
field can be thought of as arising from the tunnelling of particles between the two Rindler sectors. 



VI. CONCLUSIONS 



In conclusion, we see that the tunnelling picture is not a valid picture in general. In the case of the electric field, it 
can successfully be applied to both the time dependent and time independent gauges but such a picture is valid only to 
a limited extent. When applied naively to black hole like spacetimes it gives a temperature that is twice the expected 
value. However, by taking into account the presence of the two disjoint Schwarzchild sectors in the full Kruskal 
manifold this discrepancy can be corrected. The semi-classical prescription given in the paper takes into account the 
asymmetry between the two sides of the horizon and gives the correct result without requiring the Kruskal extension. 
An interesting aspect in the reduction to the effective Schrodinger problem is that the semi-classical analysis after the 
reduction does not give a thermal spectrum for all energies. Only when the mass of the black hole is much greater 
than the Planck mass does the spectrum reduce to a thermal form. 

But the tunnelling picture does seem to give a nice interpretation of particle production in an electric field as arising 
due to tunnelling between the two disjoint sectors of the Rindler spacetime. Though we have only given a heuristic 
argument in this paper, we will explore this issue further in a future publication. 
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APPENDIX: GENERALIZATION TO (3+l)~DIMENSIONS 



The g enera lization to (3+1) dimensions is straightforward. In this section, the semi-classical analysis given in 
section (IV A) will be briefly outlined. Further, it will be shown that the effective Schrodinger equation in (3+1) 
dimensions is the same as Eqn (|9l|). 

Consider the metrics ( p9| ) and (B0|). We will work with ( |59| ) which is in spherical polar co-ordinates. The results 
obtained are extendable to ( |60|) in a straightforward manner. 

The Klein-Gordon equation, written using the metric (59), is 



1 <9 2 $ Id 



B(r) dt 2 r 2 dr 



<9$\ 



i 



d 



dr ) r 2 sin(6>) 89 
1 



d$ 

sm(9) — 



<9 2 $ 



r 2 sin 2 (6>) dcf> 2 

Since the problem is a spherically symmetric one, one can put $ = s $>(r,t)Y l 



to obtain, 



1 d 2 ^ Id 



B(r) dt 2 r 2 dr 



dr 



1(1 + 1) m 2 



-3- * = 



(127) 



(128) 



Making the ansatz W = exp((i/h)S(r,t)) and substituting into the above equation, we obtain, 



1 



dSY 



B(r) \dt J 
h 



dS 



1(1 + l)ft 2 



B(r) dt 2 1 ; 9r 2 



1 d{r 2 B) dS 



= 



'■ dr dr 

Expanding S in a power series as in Eqn (|69|), we obtain, to the zeroth order in h/i, 

dr 



So 



-Et± 



B(r) 



yjE 2 -B(r){m 2 +L 2 /r 2 ) 



(129) 



(130) 



where L 2 = 1(1 + l)7i 2 is the angular momentum. It is easy to see from the above equation that near the horizon, 
the pres ence of the L 2 term can be neglected since it is multiplied by B(r). Therefore, the semiclassical result of 
section fl!VA|) follows even in the case of (3+l)-dimensions. The semi-classical ansatz is valid in this case as can be 
seen by calculating explictly the higher order terms in the expansion for S. All these terms have a singularity at the 
horizon of the same order as that of So and they contribute to the semiclassical propagator either as phase factors or 
as terms multiplied by powers of h which are entirely negligible. Expanding the Klien-Gordon equation for $ using 
the metric ( |60| ) gives analogous results and will not be explictly given. 

Now, consider the reduction of Eqn (127) to an effective Schrodinger equation. Setting 



* = exp(-iEt/K)Y l m (e > ^(r) 



and substituting into Eqn (127), we obtain, 



B(r) 



d 2 ^ 1 d(r 2 B) dH 



dr 2 



dr dr 



E 2 



h 2 B(r) 



n r z 



* = 



(131) 



(132) 



Making the substitution 



we get the result, 



d 2 Q 
dr 2 



1 

B 2 



(B>) 



>Y2 



* = 



--Q(r) 



h 2 



l 

B 



B^_ 



r 



k 2 r 2 







(133) 



(134) 



where B' and B" are the first and second derivatives of B(r) respectively. Near the horizon r = tq, using the expansion 
for B(r) given in Eqn (B3), it is easy to see that the 1/i? 2 term in the above equation dominates over the 1/B term 
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which can therefore be neglected. The resulting Schrodinger equation is the sa me as in Eqn (|9l|). Therefore, even in 
(3+1) dimensions, the semiclassical and quantum mechanical results of section ( 1VB ) are the same. It can be easily 
proved that the effective Schrodinger equation for the cartesian metric ([jo]) gives the same result. 
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